Maps That Take Lines to Circles, in Dimension 4 



V. Timorin 



Abstract. We list all analytic diffeomorphisms between an open subset of the 4-dimen- 
sional projective space and an open subset of the 4-dimensional sphere that take all line 
segments to arcs of round circles. These are the following: restrictions of the quaternionic 
Hopf fibrations and projections from a hyperplane to a sphere from some point. We prove this 
by finding the exact solutions of the corresponding system of partial differential equations. 

1 Introduction 

Let U be an open subset of the 4-dimensional real projective space MP 4 and V an open 
subset of the 4-dimensional sphere S 4 . We study diffeomorphisms / : U — > V that 
take all line segments lying in U to arcs of round circles lying in V. For the sake of 
brevity we will always say in the sequel that / takes all lines to circles. The purpose 
of this article is to give the complete list of such analytic diffeomorphisms. 

Remark. Given a diffeomorphism / : U — > V that takes lines to circles, we can 
compose it with a projective transformation in the preimage (which takes lines to 
lines) and a conformal transformation in the image (which takes circles to circles). 
The result will be another diffeomorphism taking lines to circles. 

Example 1. For example, suppose that S 4 is embedded in M 5 as a Euclidean sphere 
and take an arbitrary hyperplane and an arbitrary point in ]R 5 . Obviously, the pro- 
jection of the hyperplane to S 4 form the given point takes all lines to circles. The 
restriction of this projection to some open subset of the hyperplane is a diffeomor- 
phism between this subset and its image. Diffeomorphisms obtained in this way will 
be called classical projections. Of course, classical projections live in all dimensions, 
not only 4. 

Example 2. Another example comes from the quaternionic Hopf fibration. Let us 
recall the definition. Consider the standard projection from the left (resp., right) 
quaternionic 2-dimensional vector space H 2 to the left (resp., right) quaternionic pro- 
jective line HP 1 . The latter is identified with S 4 , a conformal quaternionic coordinate 
being the ratio of the homogeneous coordinates (in a specified order). Clearly, this 
projection descends to MP 7 — the real projectivization of H 2 . Thus we obtain a map 
from RP 7 to S 4 . This is the quaternionic Hopf fibration. It is known (see e.g. p^) 
that the quaternionic Hopf fibration takes all lines to circles. Therefore, the same is 
true for the restriction of it to any 4-dimensional projective subspace of MP 7 . 



Our main result is the following 

Theorem 1 Suppose that an analytic diffeomorphism f between an open set in RP 4 
and an open set in S takes all lines to circles. Then it is either a restriction of a 
classical projection or a restriction of a (left or right) quaternionic Hopf fibration. 

Remark. When we say that a map / : U C RP 4 — > S 4 is a restriction of a 
classical projection, we mean the following. There is a projective map i from RP 4 
to R 5 defined everywhere on U, a conformal identification j of 5 4 with a Euclidean 
sphere in R 5 and a central projection ir to this Euclidean sphere from some point such 
that / = j^ 1 o 7r o i. Analogously the statement "/ is a restriction of a quaternionic 
Hopf fibration" has the following meaning. There is a projective map i form RP 4 
to RP 7 , a (right or left) quaternionic Hopf fibration tt : RP 7 — * S 4 and a conformal 
identification j of the sphere in the image of / with the sphere in the image of 7r such 
that / = j 1 o tt o i. Note in particular, that Theorem^is insensitive to a projective 
transformation in the preimage and a conformal transformation in the image. 

Actually, it is not hard to see that we can always assume j to be fixed from the 
very beginning. 

History of the problem. The problem of finding maps that take lines to circles 
came from nomography (for an introduction to nomography sec 2 ). G.S. Khovanskii 
in 1970-s posed the following question: find all diffcomorphisms between open subsets 
of R 2 that take lines to circles, or, in the language of nomography, that transform 
nomograms with aligned points to circular nomograms. Circular nomograms are more 
convenient in practice whereas nomograms with aligned points are theoretically easier 
to deal with. This problem was solved by A.G. Khovanskii [3] who proved that all 
such diffeomorphisms come from projections of planes in R 3 to Euclidean spheres. 
Actually, his result was stated in a different form: up to a projective transformation 
in the preimage and a Mobius transformation in the image there are only 3 such 
diffeomorphisms, and they correspond to classical geometries. Izadi in 0] extended 
the results of Khovanskii to the 3-dimensional case. 

It turned out (see [o]) that Khovanskii's theorem does not extend to dimension 
4. The simplest counterexample is a complex projective transformation which takes 
lines to circles but does not come from a projection of a hyperplane to a sphere. This 
is a particular case of example 2 above. In it was proved that all ample enough 
rectifiable bundles of circles passing through a point in R 4 are obtained by means of 
example 2 (those obtained from example 1 can be obtained from example 2 as well). 
This may be regarded as a first step in proving Theorem^ This article completes the 
proof. 

Outline of the proof of Theorem ^ Suppose that / : U — > V is a diffeomorphism 
taking all lines to circles. First choose a projective identification of U with a region in 
R 4 and a conformal identification of V with a region in H, the skew-field of quaternions. 
Then / can be regarded as a quaternion-valued function on U. 

For any constant vector field a denote by d a the Lie derivative along a. Put A a — 
d a f. The quaternions A a are best to be understood as components of a quaternion- 
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valued differential 1-form A on U, namely, the differential of /. This form is closed: 
dA = or, in components, d a Ap — dpA a for any pair of constant vector fields a and 
/? on U. 

By Lemma 5.2 from [5], at every point of U there exists a quaternion B a that 
depends linearly on a and satisfies one of the equations 

do^A a B a A a or 3 a A a = A a B a . 

The products in the right-hand sides are in the sense of quaternions. Assume that 
the first equation holds everywhere on U (this is not very restrictive: see the end of 
Section 4). We show in Section 2 that the integrability condition for this equation is 

daBp = ^B a Bf 3 + ho(AZAfj + ~ApA a + A a ~Ap). 

where C is a quaternion-valued function. We call it the first integrability condition. 
Then we use once again that / takes lines to circles, now to deduce that C is real- 
valued. If C = identically, then the first integrability condition is integrable, and it 
gives us the maps from example 2. This will be proved in Section 3. 

Suppose now that C =/= 0. Then we need to write down an integrability condition 
for the first integrability condition — the second integrability condition. It will be 
obtained in Section 4. The second integrability condition turns out to be a differential 
relation between A, B and C only, with no additional parameters. It expresses all first 
logarithmic derivatives of C in terms of the values of A and B. The fact that C is 
real imposes a restriction on possible values of A and B. Namely, B a = p(A a ) + A a q 
where p is a real- valued 1-form and q is a quaternion. Now fix some admissible values 
of A, B and C at some point of U. We will show in Section 5 that there exists a 
map from example 1 with the same values of A, B and C at the given point. On the 
other hand, from the second integrability condition it follows that such analytic map 
is unique. Thus our map / is from example 1. 



2 First integrability condition 

In this section, we obtain an integrability condition for the equation 

d a A a — B a A a . (1) 

Since A a are derivatives, we have d a Ap = dpA a . From this and from equation (1) 
it follows immediately that 

d a Ap = ^{B a A + B () A a ). (2) 

An integrability condition for equation (1) is obtained from the equality dadpA^ — 
dpdaAj by expressing derivatives of A with the help of (2). If we introduce a tensor 
C given in components by 

Cq/3 = daBp — -B a Bp, 



3 



then the intcgrability condition reads as follows: 

C a pA 1 + C ai Ap = Cp a A 7 + CpyA a . (3) 
If we apply to equation (3) the transposition of indices [3 <-> 7, then we have: 

C a pAy + C a jAp — C^ a Ap + Cjf)A a . (3') 
Take the sum of equations (3) and (3'). It can be written as follows: 

(C/3 7 + C\p)A a = (2C Q7 - C Ja )Ap + (2C a/ 3 - Cf3 a )Ay. (3") 

Put (3 — 7 in equation (3). We obtain that 

(2C a f3 - Cp a )Aj3 = CppAa. (4) 

or, replacing (3 with 7, 

Replace 2C Q/ g — C^) a and 2C Q7 — C 7Q in (3") by their expressions obtained from (4) 
and (4'), respectively. We obtain 

C f3j + C 7/3 = CppA^A^A^A- 1 + C~ n A ,(A~ 1 Ap)A~ 1 . (5) 

Take a point in £/. In equation (5) written at this point, a, (3 and 7 may be 
arbitrary vectors. Fix (3 and 7, and let a vary. Then A a runs over all quaternions, 
since / is a diffcomorphism. We need the following 

Lemma 2 Let a and b be arbitrary quaternions, and x an imaginary quaternion ( i. e. 
not a real number). Suppose that the number ay + by~ x stays the same for all quater- 
nions y obtained from x by inner conjugations (i.e. y = qxq^ 1 for an arbitrary quater- 
nion q). Then b = a\x\ 2 . 

Proof. Denote by xo the real part of x. Then it is easy to see that 

2bx 



ay + by 1 = 



V ■ 



The second term in the right-hand side is independent of y. Therefore, the first term 
should be also independent of y which is possible only if the coefficient a — b/\x\ 2 
vanishes. □ 

We can apply this lemma to equation (5) where we put a = Cpp, b = C 77 and 
x = Ap Ay. If (3 is not parallel to 7, then x is imaginary. Thus by the lemma we have 

C/3/3 _ C 77 



\M 2 W 

It follows that the left-hand side is independent of [3. Hence 

C 0P = C\Atf (6) 
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where C is a quaternion which does not depend on the vector /?, but it may depend 
on a point — so C is a quaternion-valued function on U. 
Plug in (6) to (5): 

C 01 + C 7/3 = C(A A. ( + AyAp). (&) 

Comparing equations (6') (with (i changed to a and 7 changed to (3) and (4) we can 
conclude that 

C aP = \c(A^Ap + 7T p A a + A a ~Ap). (6") 
This is the first integrability condition. 

3 Examples 

Let us now take a closer look on examples 1 and 2 form Section 1. 

Example 1. Consider the left quaternionic Hopf fibration it : RP 7 — > S 4 . A point 
in the preimage is represented by a pair of quaternions (y, z) up to multiplication of 
both y and z by a common real factor. The image of this point under 7r is an element 
of the left projective quaternionic line with homogeneous coordinates y and z. We can 
take y~ 1 z as an afhne conformal quaternionic coordinate in the image. 

Suppose that / is obtained as the composition of tt with some projective embedding 
of (a part of) R 4 to MP 7 . This projective embedding sends a point x £ R 4 to a point 
in RP 7 with coordinates 

y = L{x), z = M{x), 

where L and M are some afhne maps from R 4 to M 4 . Thus in the given coordinates / 
looks as follows: 

J '. x 1 ► L{xy 1 M{x) (Ei) 

where the multiplication and the inverse are in the sense of quaternions. 

Denote by L and M the linear parts of L and M, respectively. This means that L 
(resp., M) is a linear operator such that L (resp., M) is a composition of L (resp., M) 
and a translation. In other words, L and M are differentials of L and M at any point. 

Compute A a = d a f for our map /: 

A a (x) = -L- 1 {x)L{a)L- 1 (x)M{x) + L- 1 {x)M{a). 
Differentiating once again along a we obtain: 

B a = -2L- 1 (x)L(a). 
We can see from here that d a B[ } = ^B a Bp, so in this case C = 0. 
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Case C = 0. Let us work out the case when C vanishes everywhere in U . In this 
case the first integrability condition (6") reads as follows: 

d a Bp = -B a B fi . 

Consider this equation together with equation (2). This is a system of partial differ- 
ential equations of first order which expresses all first derivatives of A and B through 
the values of A and B. Denote this system by S. Fix any point x G U and initial 
values of 1-forms A and B at x. Thus for any vector a at i we know A a and B a . 
Note that if a solution of S with given initial values exists, then it is unique. Indeed, 
all higher derivatives of A and B can be expressed through the initial values by means 
of the system 8. The existence of a solution follows from example 1. In this example, 
at the given point x linear maps A:aH A a and B : a^> B a can be arbitrary. 
We have just proved the following 

Proposition 3 If C = identically in U, then the general solution of equations (2) 
and (6") is given by (Ei). 

Example 2. Consider a Euclidean sphere S in R 5 . We can choose this sphere to be 
centered at the origin and to have radius 1. Introduce a conformal coordinate system 
on S by means of the stereographic projection j to the equatorial hyperplane with 
the center at the Noth pole. To this end we need to choose an orthogonal splitting 
R 5 =R 4 xB the latter factor being the line between the North and South poles. Thus 
any point in R 5 is represented as (y, z) where y G R 4 and z G R. The point on S 
corresponding to (y, 0) by means of the stereographic projection is 

|2 



3 [y) i,l + |y| 2 'l + |y|V 



Suppose that / is the central projection from a horizontal hyperplane z = z\ to S 
with the center (0, z ). Under this projection, a point (x, z\) corresponds to a point 
j^ 1 (y) on S if and only if 

2y(z - zi) 



(z + l) + (z -l)M 2 
We can now rescale y and choose specific values of zq and z\ so that 



i + M 2 ' 

This can be considered as an implicit equation defining / : x i— > y. Recall that / 
was defined as a classical projection. Another interpretation of / is that it estab- 
lishes a correspondence between the Klein and the Poincare models of the hyperbolic 
(Lobachevsky) geometry. 

Multiply both sides of equation (E2A) by the denominator to get y = (1 + \y\ 2 )x 
and then differentiate both sides along a: 



G 



Here (•, •) denotes the Euclidean inner product corresponding to the given quaternionic 
structure. Namely, for 2 vectors £ and r\ we have (£, rj) = Re(£r?) where in the right- 
hand side £ and r\ are considered as quaternions. To express A a it terms of y only we 
need to find an expression for (y, A a ). For that, take the inner product of both sides 
of equation (E 2 -2) with y. This leads to a linear equation for (y,A a ). Solving this 
equation, we obtain 

( i + M 2 ) 2 
i-M : 

Substitute (E 2 .3) into (E 2 .2). We now have a formula for A a in terms of y: 



(y,A a )= y J m J (y,a). (252.3) 



Differentiating equation (£^-4) along a we can see that 



B a = 2 2 



1 

Differentiate this equation once again along a. Then we find that 



2"° (1-M 2 ) 2 ' 
From equation (6) it follows that 

G 



C = 



(i-M 2 ) 2 



in this example. In particular, C may be nonzero. Note that in this example C is real 
everywhere. 



4 Second Integrability Condition 

The case C = has been already worked out in the previous section. Now assume 
that C ^ at a given point. A priori, C is a quaternion-valued function. But now we 
are going to prove the following 

Lemma 4 The function C is real-valued. 

This lemma does not follow form equation (1). We need to use directly the fact 
that / takes lines to circles. 

Proof. Take an arbitrary point xo £ U and a vector a at x . The line I passing 
through x in the direction of a, is mapped to a circle under /. Let t 1— > x = x(t) be 
the affine parameterization of / such that x(0) = xq and i(0) = a. 

First assume that B a is real at x . Then f(l) must be a line. Therefore, lm(B a ) 
must vanish everywhere on I. It follows that C is also real in this case. 
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Suppose now that lm(B a ) ^ 0. By Proposition 5.4 from [S] the center of the circle 
f(l) is 

f(x) - (ImB^i))- 1 ^^). 
Hence this expression does not depend on the choice of a; 6 I. Differentiate it by t: 

A a + (ImB^)" 1 (^lm{B 2 a )+lm(C)\A a \ 2 ^j (ImB Q )-M Q - {ImB^B^ = 0. 

It follows that Im(C) = 0. □ 

The first integrability condition is a partial differential equation on B. Let us try 
to write down the integrability condition for this equation. 

We will use the language of quaternion- valued differential forms. By definition, the 
algebra of quaternion- valued differential forms is obtained from the algebra of ordinary 
real-valued differential forms by taking the tensor product with the quaternions over 
reals. In particular, this definition says how to multiply quaternion-valued forms. As 
was already mentioned, A and B are thought of as quaternion-valued I-forms. 

From the first integrability condition it follows that 

dB= )-B AB + \c{AaA). (6"') 

Take the differential of the both parts of this equation, using relations d 2 B = and 
dA = 0: 

dCAAAA= ^(BA A A A- AaAaB). (7) 

Proposition 5 Suppose that the values of the 1-forms A and B and the function C 
are given at some point, and the map a A a is one-to-one at this point. If an analytic 
solution f of equations (2) with these initial values of A, B and C exists, then it is 
unique. 

Proof. Since A is non-degenerate, the operator of the right wedge multiplication 
by A A A is invertible. This can be verified by a simple direct computation. Therefore, 
equation (7) expresses all first derivatives of C through A and B. Consider this equa- 
tion together with equations (2) and (6"). We obtain a system of partial differential 
equations which expresses all first derivatives of A, B and C in terms of the values of 
A, B and C only. Therefore, by successive applications of this system we can express 
all higher derivatives of A, B and C at the given point through the initial values at 
this point. The proposition now follows. □ 

Note also that A, B and C determine the 3-jet of /. Namely, from equations (2) 
and (6") it follows that 

f _ ft \ . A . Bx-x a A x ^ Xo (f-B^-xo + C\A x - Xo \ 2 )A x _ Xo 

J — J\ x 0) + A. x - Xo -\ 1 h . . . 

In particular, Proposition can be reformulated as follows: 

Proposition 6 If there exists a diffeomorphism f : U C K 4 — > V C R 4 with a given 
3-jet at some point xq £ U such that the image of each line segment lying in U is an 
arc of a circle or a line segment in V , then such diffeomorphism is unique. 
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5 Admissible 3-jets 



From now on wc identify the space of preimage with the space of image and assume 
that the distinguished point is the origin, /(0) = and A a = a for all vectors a at 
0. This may be achieved by a linear change of variables in the preimage composed 
with a translation in the image. Thus if x = x a + ix 1 + jx 2 + kx 3 denotes the natural 
quaternionic coordinate, then A = dx = dx° + idx 1 + jdx 2 + kdx . 

Equation (7) imposes a restriction on B due to the fact that C is real. Let us find 
all admissible values of B at by solving (7) as a linear system on B and dC/C. We 
have A A A = dx A dx = iui + jcj 2 + ku 3 where u)\, lo 2 and u 3 are real 2-forms given 
in coordinates by 

Wi = 2(dx 1 A dx" + dx 3 A dx 2 ), 
uj 2 = 2{dx 2 A dx° + dx 1 A dx 3 ) , 
uj 3 = 2{dx 3 A dx" + dx 2 A dx 1 ). 

Denote by 7 the real 1-form 2dC/C and suppose that B = B° + iB 1 + jB 2 + kB 3 
where B°, B 1 , B 2 and B 3 are real 1-forms. Equation (7) can be now rewritten as the 
following system: 

2(B 2 A u 3 - B 3 A uj 2 ) = 7 A Wi, 

2(B 3 A uj x - B 1 Alu 3 ) = 7 Aw 2 , (8) 
2(B X A uj 2 - B 2 A uji) = 7 A w 3 . 

Denote by B^ and j v (u = 0, 1, 2, 3) the components of the 1-forms B^ and 7 so that 

3 3 
B» = Y J Byx", 7 = ^7^. 

v=0 i/=0 

System (8) yields the following equations: 

2(-B 2 + B 3 ) = - 7l , 2(-B 2 - B 3 ) = 70, 2(B 2 - B 3 ) = -73, 2(B 2 + B 3 ) = 72 , 
2(-B 3 + S3 1 ) - - 72 , 2(B 3 + fii) - 73 , 2(-B 3 3 - B\) = 70, 2(B 2 3 - ^) = - 7l , 
2(-B 2 1 + B?) = _ 73j 2 {B\ - B 2 ) = - 72 , 2(B^ + B 2 ) = 7l , 2{-B\ - B 2 2 ) = 7 o- 

Exclude 7^ from these equations: 

B l = B 2 = B 3' 

B h = B l = ~ B li 
B 2 = B\ = -Bl 
B 3 = B\ = -B\. 

The relations displayed above mean that B a = p(a) + aq where p is a real-valued 
1-form and q is a quaternion. Hence the 3-jet at of / should be as follows: 

(p(x) + xq)x (%(p{x) + xq) 2 + C\x\ 2 )x 
X+ 2 + 6 ' [) 

Proposition 7 For an arbitrary real-valued linear function p on M 4 and an arbitrary 
quaternion q, a local diffeomorphism with 3-jet (9) exists. Moreover, it can be chosen 
to be one of the classical projections. 
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Proof. First consider a map given by the implicit equation 

- y 

x i-fM 2 ' 

This is clearly a map from example 2. It has the 3-jet x + ^\x\ 2 x. This is a partial 
case of (9) where p = q = 0. To achieve other given values of p and q, compose this 
map with the following Mobius transformation: 



It is easy to see that the composition has 3-jet (9). □ 



Concluding remarks. We should now make things add up in the proof of theorem 
□ 

We have always assumed that d a A a = B a A a the multiplication by B a being 
from the left. The case of the right multiplication is completely analogous. The only 
difference is that for C — we would have the right quaternionic Hopf fibration instead 
of the left one. In general there is an open subset U' of U such that the multiplication 
by B a is either form the left everywhere on U' or from the right everywhere on U'. 
Therefore theorem holds on U'. By the uniqueness theorem for analytic functions it 
holds then on U . 

If C = everywhere on U, then by Proposition we have a quaternionic Hopf 
fibration. Otherwise C nowhere vanishes on some open subset of U. Then by Propo- 
sitions and the map / is a classical projection on this subset. By the uniqueness 
theorem / is a classical projection on U. 
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